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, Abstract 

' We are analyzing the hydrodynamics of 200 A GeV S+S colhsions using a new 

approach which tries to quantify the uncertainties arising from the specific im- 
' plementation of the hydrodynamical model. Based on a previous phenomeno- 

, logical analysis we use the global hydrodynamics model to show that the amount 

' of initial flow, or initial energy density, cannot be determined from the hadronic 

, momentum spectra. We additionally find that almost always a sizeable trans- 

' I ' verse fiow develops, which causes the system to freeze out, thereby limiting the 

O , flow velocity in itself. This freeze-out dominance in turn makes a distinction 

between a plasma and a hadron resonance gas equation of state very difficult, 
whereas a pure pion gas can easily be ruled out from present data. To complete 
^ ^ ' the picture we also analyze particle multiplicity data, which suggest that chem- 

, ical equilibrium is not reached with respect to the strange particles. However, 

\^ • the overpopulation of pions seems to be at most moderate, with a pion chemical 

5^ , potential far away from the Bose divergence. 

1 Introduction 

The success of hydrodynamics as a model for the space-time evolution of heavy- 
ion colhsions at the BEVALAC with center-of-mass energies of several hundred 
MeV/nucleon |, || suggests its basic apphcability to the ultrarelativistic heavy- 
ion colhsions at the Brookhaven AGS with ~ 5 A GeV or at the CERN SPS with 
^ 20 A GeV (e.g. |, g |, |, §). Since the hydrodynamical model |, ||, |n| 



describes the matter during the collision in terms of the collective variables energy 
density e, pressure P and flow velocity n'^, its validity should actually improve at 
higher energies because of the strongly increasing particle multiplicities. Although 
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the model's primary task is the description of the inclusive hadron spectra, which 
are measured in several experiments at the AGS |12, 13] and the SPS |13, 14], its 
reasonably precise picture of the collison dynamics can furtheron serve as a basic 
framework for the computation of almost any other observable, e.g. J/'0, strangeness, 
p, dileptons, strangelets, etc. 

Moreover, such an effort is also accompanied by the hope that once a reliable hy- 
drodynamical description is achieved we will automatically have information about 
the yet unknown equation of state of hadronic matter, which serves as the basic 
input to hydrodynamics. Especially, a first order phase transition from normal 
hadronic matter to a quark gluon plasma, as suggested by lattice QCD calculations, 
would result in a coexistence region with dP/de — > and, in principle, might lead 
to qualitatively new dynamical signatures. 

However, we encounter the problem that some of the assumptions underlying the 
hydrodynamic model may not be sufficiently accurate in ultrarelativistic collisions 
with center-of-mass energies one order of magnitude higher than at the BEVALAC. 
In particular, the question whether and when local thermalization can be reached 
is still under debate, and we do not attempt to answer it in this paper. Instead 
we will take sufficient thermalization for granted some time after the first hard 
collisions and only assume validity of the hydrodynamic model thereafter. Apart 
from thermalization, our approach will differ considerably from previous publications 
]|3|, 1^, |5|, |6|, 1^ in that we will try to resolve the other important uncertainties arising 
from the specific implementation of hydrodynamics. To achieve this goal we are 
careful to test each hypothesis when it is introduced, leading to the organization of 
the paper as follows. 

Based on a phenomenological analysis of the hadronic spectra from S+S at 
200 ^GeV [H], we extract in Sec. |2| from the measured particle spectra as much in- 
formation as possible about the hydrodynamical behaviour at freeze-out. In Sec. ^ 
we use global hydrodynamics as a computationally efficient method ]|l^ , to investi- 
gate to what extent the initial hydrodynamical state is constrained by the data and 
find that the amount of initial longitudinal flow cannot be decided from the data 
alone. A closer look at two possible extreme scenarios in Sec. ^, namely Landau- 
and Bj or ken-like initial conditions, shows striking similarities in their temperature 
and transverse flow at freeze-out. In Sec. |5| we show that the transverse expansion 
is the most important contributor to the freeze-out process and is thus invariably 
limited by itself. In Sec. ^ we confront our model with the experimental data and 
note good agreement for a variety of underlying assumptions regarding entropy gen- 
eration and freeze-out description. Additionally, we test different equations of state 
and find that only the pure pion gas can be ruled out while the distinction between 
the plasma and the hadron resonance gas by hydrodynamical means is very hard 
to accomplish. In Sec. ^ we finally address the issue of chemical equilibration and 
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the pion chemical potential by computing particle ratios and multiplicities. Sec. ^ 
summarizes our results and conclusions. 

2 The data from a hydrodynamical point of view 

The measured single particle distributions confer basically two distinct types of 
information: The shape of the measured spectra reflects the temperature and flow- 
dynamics of the system, while the absolute normalization of the spectra is connected 
additionally with the size of the collision zone and the amount of chemical equili- 
bration among the particle species. 

For the moment we will focus exclusively on the shape, leaving the normalization 
arbitrary until Sec. |^ where we return to this issue in detail. Since at SPS energies 
pions are the overwhelming majority of all produced secondaries, with kaons, nu- 
cleons, strange baryons and antibaryons following in descending order, we will base 
our analysis preferentially on the pions. This ensures on the experimental side that 
the spectra can be measured quite accurately with small statistical errors while on 
the theoretical side the understanding of different effects on the spectra has been 
developed the furthest. In any case, because the other particles seem to follow the 
pions (Fig. |2|), all species would give very similar results. 

As shown in the dynamics of the collision zone shows up in the hadronic 
spectra, as measured by the NA35 collaboration |17, in two distinct ways: 



The longitudinal flow determines the rapidity distributions, while the transverse flow 
influences the transverse mass spectra. 

The longitudinal flow can be extracted most accurately from the pion rapidity 
spectrum 0, 120, |l5| . The extracted value is supported by the spectra of all other par- 



ticles measured by NA35 |17|, y^, |19[. The rapidity spectra of all produced particles 
are well described by a boost invariant distribution around the center of the collision 
zone at ycm = 3, with the flow rapidities limited to the interval (— //max; +^max) with 
r?^ax = 1.7±0.2 m (see Fig. IT]). 

The situation is less clear for the transverse mass spectra, which are affected both 



by resonance decays |21| and transverse flow. We found that for any given radial 
fluid velocity < /?s < 0.7 c it is always possible to find a temperature T{l3s), which 
fits the spectra for all particle species simultaneously [^] (see Fig. |2|). The reason 
is that with transverse flow the spectra can in good approximation be described in 
terms of an apparent or blueshifted temperature [^] 

(1) 

from which the temperature T and average radial flow (Pr) cannot be easily sepa- 
rated. 
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Fig. 1: Every point (T, 77max) on these 
lines corresponds to a fit of the com- 
puted y-spectra (including resonance 
decays and longitudinal flow) to the 
measured ones The longitu- 

dinal flow T^max ~ 1.7 ±0.2 fits all parti- 
cle spectra almost independently of the 
temperature T, with exception of the 
protons, which still carry a big amount 
of their initial longitudinal motion . 




Fig. 2: Every point (T, /3s) on these 
lines corresponds to a fit of the com- 
puted TO^-spectra (including resonance 
decays and transverse flow) to the mea- 
sured ones 0, Temperature T 
and transverse flow f3s have similar ef- 
fects on the spectra so that /3s can not 
be extracted unambiguously ITs). 



The amount of transverse flow thus can not be determined from the particle 
spectra alone. We have shown [^] that the method of global hydrodynamics fl^ 
combined with a reahstic freeze-out criterium is well suited to give "circumstantial 
evidence" for a transverse flow of Ps ~ 0.5c |23|. Here we will present a more 
thorough theoretical analysis to gain a better understanding of the dynamics of the 
collision zone. The method of global hydrodynamics is ideal for such an analysis, 
because on the one side it contains all necessary elements for a consistent description 
of the dynamics (coupling of longitudinal and transverse flow by the hydrodynamic 
equations as well as a dynamical freeze-out criterium) and on the other side, due 
to its numerical effectiveness, it allows for an in-depth phenomenological analysis 
including extensive parameter studies. 



4 



3 Extracting the Initial Conditions 



Hydrodynamics connects the thermalized initial state with the frozen- out final state. 
With its help the freeze-out conditions (Tab. |^), reflected in the observed particle 
spectra, can be related to the initial conditions and provide new insights into the 
dynamical evolution. Some of the initial conditions are known from the size of the 
projectile nucleus and the measured energy loss of the projectile protons [24, 18]. 
At the freeze-out point we know the maximal fluid rapidity (Fig. |^) as extracted 
from the y-spectra. The m^-spectra only yield a correlation between the transverse 
fluid velocity (3s and the local temperature T at freeze-out (Fig. |2|), which is rather 
difficult to exploit and the discussion of which we postpone until Sec. ^. In principle, 
it is also possible to determine the final radius using pion interferometry (HBT) p5[| . 
Unfortunately, up to now the theoretical and experimental uncertainties in the ana- 
lysis of the 2-pion (and 2-kaon) correlation functions are still too large (> 0.5 fm) 
to render them very useful for our purpose [26|. 



quantity 


initial state 


from: 


final state 


from: 


R 


4fm 


projectile size 


? 


(HBT) 







defined 


(A,T) 


dn 
mTdiriT 


Z 


from Etot 


^loss = 313GeV 


? 


? 


V 


? 


? 


1.7 


dn 
dy 


£ 


? 


? 


? 


? 



Tab. 1: The initial and final state parameters of the global hydrodynamics are 
only partially known. Through hydrodynamics we can combine our knowledge 
about the initial and final state and constrain the unknown parameters. 



The idea of this paper is to use the global hydrodynamics framework of Ref. |16] 
to fill in (most of) the question marks in Tab. |l|by establishing a connection between 
the initial and the final state. We proceed strategically and first narrow down the 
uncertainties in the initial conditions by choosing from all combinations of initial 
longitudinal flow ryo and initial energy density eq only those pairs (%)£o) which 
after hydro dynamical evolution produce at freeze-out a longitudinal flow value of 
r]f = 1.7. In Fig. ^ we show the result for a hadron resonance gas equation of state 
with speciflc entropy S/A = 30. 
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For the case of complete stopping (no longitudinal flow at the beginning, r]Q = 0) 
we need very high energy densities eo > 10 GeV/fm^, so that the strong longitudinal 
pressure gradient can accelerate the matter to rjf = 1.7 at freeze-out. For larger 
initial flows considerably smaller initial energy densities suffice to reach the same 
freeze-out value. In the extreme case of very large initial flows r/o ~ ^/ the system 
is already close to freeze-out initially, i.e. the initial energy densities are very low, 
and the hydrodynamics barely generates additional longitudinal flow. 
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Fig. 3: Shown is the initial en- 
ergy density eo(''7o) which, for a given 
initial longitudinal flow 770, generates 
at freeze-out the measured flow 77/ = 
1.7. The hydrodynamics is based on a 
hadron resonance gas EOS with S/A — 
30. 



Fig. 4: Temperature T and transverse 
flow Psj at freeze-out as determined 
from the longitudinal flow 77/ = 1.7 via 
the global hydrodynamic model. Both 
are almost constant over a whole range 
of initial conditions for the longitudinal 
flow, rjo < 1.6. 



Knowing the possible pairs of initial conditions (eO)%)) we can now study 
the systematics of the hydro dynamically generated transverse flow at freeze-out, 
parametrized by /?<j, within the allowed range of longitudinal flow parameters. Sur- 
prisingly, for a large region of initial conditions, ranging from complete stopping to 
considerable initial transparency (0 < r/o < 1-6), more or less the same (sizeable) 
amount of transverse flow fig ~ 0.5c is generated [23|. Above ryo ^ 1-6 no real 
hydrodynamical evolution takes place since the matter is already close to freeze-out 
and is accelerated very little. Although strictly speaking we cannot exclude this 
possibility within the framework presented here, we will disregard this region from 
now on, because there the assumption of local equilibrium which underlies the hy- 



drodynamic model breaks down, and a non-equilibrium description |27] appears to 
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be more suitable. 

The temperature at freeze-out turns out to be similarly insensitive to the initial 
longitudinal flow. This is largely a consequence of the constancy of Psj, since the 
temperature at freeze-out is coupled to the transverse velocity (see Sec. via the 
scattering time scale Tsca{T) and the expansion time scale Te^piR, Z,r]f , I3sj) such 
that higher velocities cause the system to freeze out already at higher temperatures 
and densities. 

From this we can already conclude that with hadronic observables alone the 
initial energy density and the initial longitudinal flow can not be determined. All 
measurable quantities at freeze-out are essentially independent of r/o (and thus of Eq). 
Other obervables, for example lepton pairs |28|, might prove useful to determine the 
initial energy density: due to their small cross sections they freeze out immediately 
and thus the measured spectrum receives contributions from the whole space-time 
region of the collision zone, with the hot early stages dominating the spectrum due 
to the T^-dependence of the Stefan-Boltzmann law for thermal radiation. 



4 Full stopping vs. partial transparency 

We will postpone the application of these results for the moment in order to analyse 
in this section in more detail their dynamical origin. We pick two extreme cases 
of the initial longitudinal flow which we call, according to their resemblance to the 
Landau and the Bjorken hydrodynamics, the L-scenario {rjQ = 0) and the B-scenario 
(770 = 1-3), respectively. Their initial and freeze-out values are given in Tab. |2|. 





L-scenario 


B-scenario 


initial 


final 


initial 


final 


R (fm) 


4 


4.71 


4 


4.57 


Ps (C) 





0.58 





0.50 


Z (fm) 


0.23 


4.0 


2.3 


7.2 


V 





1.7 


1.3 


1.7 


e (GeV/fm^) 


13.8 


0.36 


1.4 


0.26 


T (MeV) 


313 


167 


209 


157 


S/A 


30 


30 


30 


30 



Tab. 2: Initial and freeze-out values of the L and the B-sccnario, which 
are distinguished by their initial longitudinal flow 770, as computed with 
our model from the known initial values (Tab. 0) and the choice S/A = 30 
for a hadronic EOS. 



The L-scenario generates the largest transverse fiow, as one would have intu- 
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itively expected. But the excess of transverse flow compared to the B-scenario is 
rather small; at first sight this is quite puzzling, considering the large difference in 
the thermally generated longitudinal flow rjf — tjo, which is 1.7 for the L-system and 
only 0.4 for the B-system. 

In Fig. |5| we compare the dynamical evolution of both scenarios until freeze-out 
(where the lines end). The initial time of the B-scenario was displaced by hand in 
such a way that at that point the temperatures in the two scenarios agree. One 
realizes that the production of longitudinal flow is strongest at the beginning and 
then levels off considerably, even when taking into account the exponential character 
of ij. 
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Fig. 5: Comparison between 
the L- and B-scenarios shows 
remarkable similarities if one 
shifts the initial time of the B- 
dynamics so that its tempera- 
ture agrees with the L-scenario 
at some later time. The differ- 
ence in build-up of transverse 
flow is compensated by a later 
freeze-out in the B-scenario. 



t (fm/c) 



Neglecting small deviations between both scenarios - which incidentally might 
also be caused by the slightly differing longitudinally comoving coordinate systems 
in both cases - the main difference clearly lies in /J^. The L-scenario generates trans- 
verse flow from the beginning, so that there is already a flow velocity of (3s ~ 0.3 c 
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when the B-system starts. In the further evolution the L-system freezes out earher 
whereas the B-system nearly catches up in transverse flow during the remaining 
part of its expansion history. The reason for the similar transverse flow velocities 
at freeze-out thus appears to lie more in the freeze-out criterium than in the initial 
conditions. 



5 Longitudinal and Transverse Freeze-Out Time Scales 

We investigate this question in more detail in Fig. |6| by analyzing the dynamically 
determined expansion time scale Texp- Initially the evolution is strongly dominated 
by longitudinal expansion, resulting in small values of Tgxp. However, the length 
of the system quickly increases, reducing the longitudinal velocity gradients and 
thus weakening the longitudinal shear forces which could break the local thermal 
equilibrium. A simple estimate in the limit of boost invariant longitudinal expansion 
gives an expansion time scale Texp = t which grows linearly with proper time. 




1 2 3 4 5 
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Fig. 6: The expansion time scale (here 
for the B-scenario with 770 = 1-3) is ini- 
tially dominated by longitudinal expan- 
sion. After a few fm/c the radial ex- 
pansion takes over. The net result is 
a rather weak time dependence of the 
total expansion time scale. 



2 

■sea 



5 
2 







71 




S/A = 20 






S/A = 30 






- - - S/A = 40 





100 125 150 175 



T (MeV) 

Fig. 7: The scattering time scale of 
pions in a hadron gas varies with S/A 
because of the different baryon con- 
tent. Its exponential decrease at higher 
temperatures, combined with the near 
constancy of Texp, results in freeze-out 
temperatures which are always close to 
150 MeV. 



The transverse radius, on the other hand, stays for a long time close to its 
initial value. Due to quadratic i?-dependence of the transverse area of the cylin- 
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der, already rather small transverse velocities of order 13s ~ 0.4 c cause appre- 
ciable dilution and lead to transverse expansion time scales which can compete 
with the longitudinal one. This is illustrated by the non-relativistic estimate 
''"cxp ~ (n + \)(3s/R ~ 3fm/c, with n = 2 for a quadratic transverse velocity profile 
and R = 4fm/c for a sulphur nucleus. We conclude that the rarefaction at later 
stages is strongly dominated by the transverse flow, which thereby limits its own 
growth. A comparison with the scattering time scale Tsca in Fig- 0; which constitutes 



our freeze-out criterium Texp = ''"sea [l29| , 16 1, then shows why always nearly the same 
freeze-out temperature Tf is generated. 

Our whole analysis is based on this dynamical freeze-out criterium, which explic- 
itly models the competition between the local expansion, which disturbs equilibrium, 
and the scattering processes, which restore equilibrium. If we had implemented the 
freeze-out through a static criterium, e.g. via a freeze-out temperature Tf = const, 
(typically 140 MeV), this kind of analysis would not have been possible. It is inter- 
esting to observe, however, that the result of dynamical freeze-out nearly coincides 
with such a criterium, due to the rapid cooling of the system by transverse expansion 
and the exponential temperature dependence of the scattering time scale. 

In summary we arrive at the remarkable result that there is no strong correlation 
between produced longitudinal and produced transverse flow, because the former is 
mainly generated initially while the latter results from the late dynamical stages be- 
fore freeze-out. This also explains retrospectively why, in spite of the big differences 
with regard to the longitudinal dynamics, the spherical model of Ref. (3^ yields 
similar values for the transverse flow at freeze-out. 



6 Comparisons with the Data 

We can now confront our theoretical prediction for the relationship between (3s and T 
at freeze-out with the pairs (Tf,f3sj) we have extracted from the phenomenology of 
the transverse mass spectra []l5| . The hope is that this comparison permits to narrow 
down the amiguity between thermal and collective motion from the phenomenology 
of the transverse momentum data [|l^]. However, even though the phenomenological 
analysis of the data gives an anticorrelation between Tf and Psj according to Eq. (|T|), 
whereas the freeze-out criterium correlates T and Ps,f directly (faster expansion 
causes freeze-out at higher T), it is a priori by no means guaranteed that the two 
approaches yield comparable (T,Ps) values at all: the only input for the theoretical 
hydrodynamical model are the rapidity distributions, while its results are checked 
by using the independent data from the transverse mass distributions. 

For the comparison we combine the theoretically determined Psjivo) ^^i^ '^fivo) 
into a line {Tf,f3sj), which we parametrize by the unknown initial longitudinal 
flow rjo, covering the range from to 1.6 as indicated next to the line in Fig. ^. 
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Fig. 8: A comparison of 
the theoretical correlation be- 
tween temperature and trans- 
verse flow at freeze-out with 
the experimental correlation 
deduced from the pion spectra 
(Fig. |2|). The intersection at 
Ps ~ 0.5 c can be interpreted as 
evidence for transverse flow in 
the experiment. The hydrody- 
namic system with w 10% en- 
tropy generation (dotted line, 
|lq|) gives a similar result. 



From Fig. |^ we already know that the pairs {Tf,(3sj) for all initial conditions with 
genuine hydrodynamical evolution, i.e. for initial flow rapidities < ?7o < 1.6, are 
concentrated in a small region near Ps ~ 0.5 c and T ~ 150 MeV. 

It is gratifying to see in Fig. ^ that the (T, Ps )-curve extracted from the data 
crosses this region; this implies that the phenomenological (T, /3s)-correlation is con- 
sistent with hydrodynamical evolution and freeze-out, and that our dynamical pic- 
ture of the collision region successfully explains the data. We can use this agreement 
between theory and data as support for our assumption of local thermalization and 
hydrodynamical evolution and interpret the result as a theory-based proof for the 
existence of transverse flow at freeze-out. To assess the credibility of such an inter- 
pretation and to venture into new grounds, we now proceed to examine our analysis 
under various additional aspects. 

The interest for transverse flow in heavy-ion collisions has been fueled also by 
the hope to find an indicator for a phase transition into the quark gluon plasma 
1 31]. In an intuitive scenario the initial collision energy in longitudinal direction 
should be partially converted into the latent heat associated with a first order phase 
transition. The subsequent cooling and decay of the plasma would then release the 
energy isotropically to all particles, so that in particular the transverse momenta 
would profit. 

The possible existence of a phase transition to a quark-gluon plasma state in- 
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Fig. 9: Sensitivity of the 
transverse flow to different as- 
sumptions: The reference curve 
from Fig. || is labelled by S/A = 
30. Other hadron resonance 
equations of state with S/A = 
20 and 40 result in almost the 
same transverse flow /J^ ~ 0.5 c, 
whereas freeze-out at constant 
temperature T = 140 MeV 
leads to a slightly larger trans- 
verse flow. 



(c) 



fluences the hydro dynamical evolution through the corresponding modification of 
the equation of state of the expanding matter. We have tested with the global hy- 



drodynamics quantitatively the equations of state depicted in Fig. 10. They are 



presented in terms of energy density e vs. pressure P, the variables entering the 
hydrodynamical tensor T^^^ . We repeat our procedure to obtain (T,Ps) correlations 
leaving the unknown initial r/o open. 

We begin by discussing the hadron resonance gas EOS, which describes a mixture 
of relativistic ideal gases of a multitude of hadron species (tt, rj, p, w, (/>, , , p, n, 
A, A, S, H, Q, and their antiparticles) . Previously we fixed the entropy per baryon 
at S/A = 30, to achieve agreement with the measured p/n' ratio (Sec. ^). We now 
also include the choices S/A = 20 and 40, corresponding to systems with larger and 
smaller baryon densities, respectively. Clearly, all three EOS are very similar in -P(e) 
such that the baryon content has little influence on the dynamics. However, due to 
the large pion-nucleon cross section the baryon content does influence the scattering 
time scale (Fig. 0). The resulting variations in the freeze-out point can be seen in 
Fig. |. 

We have also implemented the simple freeze-out criterium Tf = 140 MeV = 
const., as it is often used in hydrodynamics calculations. This also permits us to 
compare below the resonance gas and plasma equations of state, since for the latter 
the computation of a scattering time scale for the final state hadrons is not really 
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Fig. 10: The equation of state in- 
fluences the hydrodynamics by the de- 
pendence of the pressure P{e) on the 
energy density. There are significant 
differences between a soft hadron res- 
onance gas (sohd), the plasma equa- 
tion of state including a phase transi- 
tion (dotted) and the hard pure pion 
gas (dashed). 
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Fig. 11: The transverse flow from 
a plasma EOS differs only little from 
that of a hadron resonance gas when the 
same freeze-out criterium is employed 
(Fig. H). For a pure pion gas EOS all 
initial conditions with genuine hydro- 
dynamical evolution (770 < 1.5) gener- 
ate too much transverse flow and fail to 
reproduce the data. 



meaningful except near the end of the collision. Because the value of 140 MeV 
is below the dynamical freeze-out temperatures of our standard scenario, the flow 
values in this case tend to be somewhat higher. 

An extreme, but popular EOS is that of a gas of non-interacting pions. It is 
much harder (Pnie) 3> fhad(e)) than the other equations of state and very close 
to ideal massless boson limit P = ^e. It thus generates a very large transverse 
flow (Fig. |l^), resulting in a known (e.g. ^]) disagreement with the observed rriT- 
spectra. This EOS can thus be excluded by our analysis for all hydrodynamic initial 
conditions. 

Finally we study a plasma equation of state with a first order phase transition, 
where during the mixed phase plasma and hadron gas coexist. As a representative 
from this class we have taken the EOS of U. Ornik |32], which has been obtained 
by interpolation of data points from lattice gauge theory calculations. The plasma 
EOS exhibits a bend in the coexistence region where the velocity of sound Cg'^ = ^ 
approaches zero. The comparison of Figs. ^ and 11 for Tf = 140 MeV shows that 
the plasma EOS generates a similar amount of flow as the hadron EOS. This results 
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from the fact that the transverse flow is generated at later stages, when the system 
is no longer in the plasma state but already in the hadron gas phase. 



7 Particle Ratios and Multiplicities 

While the thermal description of the momentum spectra follows from the assump- 
tion of local thermal equilibrium, which requires only elastic scattering, a thermal 
description for the relative particle multiplicities (e.g. as in |3^) implies also chemi- 
cal equilibration which requires additionally inelastic scatterings. While the former 
seems to be reasonably valid for heavy- ion collisions, the latter is much harder to 



reach in purely hadronic systems, especially for strange particles [34|. In fact, our 
approach does yield characteristic discrepancies between theoretical and experimen- 
tal particle ratios, which indicate that full chemical equilibrium has not been reached 
in 200 ^GeV S-hS collisions. 

Since particle ratios in a grand canonical description are parametrized by chem- 
ical potentials, we now introduce and investigate the dependence on the chemical 
potentials fib for the baryon number and fig for the strangeness. While fig can be 
fixed by the condition of local strangeness neutrality {ps = 0, corresponding to the 
initial state of the collision), fi^ at freeze-out can be adjusted by fitting the ex- 
perimental p/tt ratio. We prefer S/A instead of Hb as a parameter, because it is an 
invariant of the hydrodynamical evolution and thus can characterize the EOS during 
the whole expansion. 



experiment 


global hydrodynamics 




NA35 
total y 




S/A = 20 
Tf « 150 MeV 


S/A = 30 
Tf w 162 MeV 


S/A = 40 
Tf 166 MeV 


p/h 


0.28 ± 0.04 


p/tt 


0.45 


0.27 


0.21 


K^/h- 


0.10 ±0.02 




0.19 


0.20 


0.21 


A/h- 


0.080 ±0.010 


A/vr- 


0.17 


0.13 


0.11 


A/h- 


0.015 ±0.005 


A/tt- 


w 0.005 


w 0.015 


w 0.024 



Tab. 3: The particle ratios from our model in comparison with the data from NA35 ]l7| , 
|l9|. The actual number of tt~ in the experiment could be smaller than the measured 
negative hadrons by ~ 5% because of and p contributions, whereas the protons 
are measured as the difference between positive and negative hadrons. 
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The particle ratios (including resonance decays) from global hydrodynamics are 
shown in Tab. ^ in comparison with data from NA35. We compare the total mul- 
tiplicities, thereby averaging over the rapidity distributions, which may - especially 
in case of the protons - be different between our model and the experiment (Fig. 5 
in |jl5[). Again the hydrodynamics results are almost independent of the initial lon- 
gitudinal fluid rapidity r/o- As expected, the p/vr" ratio varies strongly with S/A, 
and the measured ratio can be reproduced quite nicely by the choice S/A = 30. If 
we would try to account for the reduced baryon content at central rapidities, which 
is, as mentioned, not fully reproduced by our model for the central zone, we would 
opt for a values of S/A = 40. 

The relative multiplicities of the strange particles are generally overestimated 
by the model, which at first sight seems to indicate that in the experiment chem- 
ical equilibrium for strangeness has not been established yet by the ss generating 
processes. In all three cases the value of K^/ir" is almost the same, because the 
higher kaon abundance at higher temperatures is offset in the denominator by the 
pion production from increased resonance decay contributions. The A depend very 
strongly on the freeze-out temperature and the baryon content. 

However, the apparent overprediction of the strange particle ratios relative to the 
negatively charged particles could also be due to an underprediction of the total pion 
yield by the thermal equilibrium model. As a number of authors have pointed out 
j35[ |, the pion creation processes might dominate the reverse reaction for some time 
leading to an overabundance of pions compared to their chemical equilibrium value. 
Such non-equilibrium features can be modelled by introducing a finite pion chemical 
potential It would attain equal (positive) values for particles and antiparticles, 
in contrast to the vanishing chemical equilibrium value for bosons without conserved 
quantum numbers. 

Tab. ^ shows that global hydrodynamics lacks indeed both in the number of 
pions and their energy. This lack of pion multiplicity accounts for a large frac- 
tion of the strange/non-strange discrepancy between theory and experiment seen in 
Tab. 1^. Therefore the strange sector in S+S collisions is actually much closer to 
chemical equilibrium than first suspected [^]. This was recently confirmed by a 
more thorough chemical analysis of the data in [p^]. 

By computing the total energy per pion E-j^/N.,^, we deduce from the relatively 
good agreement and the high value of Stot/^ir that the discrepancy is not so much 
a result of the shape of the pion spectrum but merely of the general underrepresen- 
tation of the pions in our model relative to strange particles. The observed small 
deviation in the ratio Et^/N-,^ could be easily due to the pions in the lowest bin 
{rriT — ni-n^ < 100 MeV/c^), which strongly influence the multiplicity because of the 
exponential character of the m^-spectrum, and which we cannot reproduce in full 
detail (Figs. 2&7 in |l|). 
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experiment 


global hydrodynamics 


NA35 


S/A = 20 


S/A = 30 


S/A = 40 




93 ±5 


55 ... 59 


62 ... 67 


64 ... 71 


E^- (GeV) 


71 ±7 


45 


54 


57 


E^-/N^- (GeV) 


0.76 


0.76... 0.84 


0.81... 0.87 


0.80... 0.89 






7.4 


6.5 


6.2 



Tab. 4: The total multiplicity and energy of tt~ are not very well represented by our 
model. However, the energy per pion is described rather well. The total entropy per 
pion is considerably higher than the value of ~ 4 for a pure pion gas, reflecting the large 
abundance of heavy particles in the resonance gas. The data are from M, M. 



If we were to consider that all the missing pions are generated by strong Bose 
condensation, we would find (using Boltzmann distributions only) a pion chemical 
potential of /x^r = 50 ... 80 MeV as a rough estimate, without taking into account 
the effect of a correspondingly harder EOS on the hydrodynamics. This value of 
/^TT is considerably smaller than = 126 MeV in [^5[ and thus sufficiently far away 
from the Bose divergence at = 139 MeV so that it does not lead to a substantial 
overpopulation at low momenta. 

8 Summary and Conclusions 

We started from a previous analysis [15| of hadronic transverse mass and rapidity 



distributions from 200^ GeV S+S collisions to estimate the maximum information 
content a hydrodynamic analysis can possibly deliver. Using all available information 
on both the initial conditions as well as the freeze-out conditions, we still cannot 
determine the initial longitudinal flow, or equivalently, the initial energy density. 
However, we find that regardless of the initial flow all possible initial conditions 
finally lead to the same temperature and transverse flow at freeze-out, which we 
model with a dynamical freeze-out criterium. This is apparently a consequence of 
the dynamical freeze-out in conjunction with the transverse expansion, which in the 
late stages of hydrodynamics proves to be much more important for the rarefaction 
than the longitudinal expansion. 

A comparison with our phenomenological flow analysis shows then good agree- 
ment with the data for a wide variety of initial conditions and freeze-out criteria. 
However, while both hadron resonance gas and plasma equation of state produce 
good agreement with the data, a purely pionic equation of state generates too much 
transverse flow and can thus be excluded for all hydrodynamic initial conditions. 
Furthermore, from the fact that the data show a lower number of strange particles 
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than our model with chemical equilibrium would produce, we deduce that chemi- 
cal equilibration has not been reached. On the other hand we determine that the 
chemical non-equilibrium of pions is not as big as to allow for a close to the pion 
mass, wc estimate instead //jr ~ 50 . . . 80MeV. 
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